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Abstract. J/ijj production at forward and midrapidity at the Relativistic Heavy-Ion Collider (RHIC) is 
calculated within a previously constructed rate-equation approach accounting for both direct production 
and regeneration from c and c. The results are compared to experimental data. The observed stronger 
suppression at forward rapidity can be qualitatively explained by a smaller statistical regeneration com- 
ponent together with stronger cold nuclear matter induced suppression compared to midrapidity. The Xc 
over J/ip ratio and ip' over J/ip ratio are also calculated. 
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1 Introduction 

The suppression of J/ip mesons in ultrarelativistic heavy- 
ion collisions, as a result of color Debye screening, has 
been suggested as a signature [I] of Quark-Gluon Plasma 
(QGP) formation a long time ago and is indeed observed 
in both Pb-Pb collisions at the CERN Super Proton Syn- 
chrotron (SPS) [2] and in Au-Au collisions at BNL Rela- 
tivistic Heavy-Ion Collider (RHIC) |3J. However, it is puz- 
zling that the observed suppression is very comparable at 
SPS and RHIC energies, since the energy density of the 
medium at RHIC is expected to be much higher than at 
SPS. The statistical model 4J, or kinetic approaches 
[B], explain this puzzle by considering the regeneration of 
charmonium from c and c quarks: at RHIC energies more c 
and c quarks (relative to charmonium) are produced than 
at SPS energies and therefore the regeneration of charmo- 
nium at RHIC largely compensates the expected stronger 
suppression (which was, in fact, predicted in Ref. [7]). Re- 
cent RHIC data [3] suggest another puzzle: charmonium 
suppression observed at forward rapidity (\y\ €[1.2,2.2]) 
is stronger than at midrapidity (\y\ <0.35), despite the 
energy density of the medium at forward rapidity being 
presumably smaller which should lead to weaker suppres- 
sion. 

J/ip production at forward rapidity has been investi- 
gated by several theoretical models. In the kinetic recom- 
bination model [5] the rapidity and transverse momentum 
distributions of J/ip are obtained through solving a rate 
equation with the gain term accounting for the continu- 
ous formation process of the J/ip's from c and c quarks 
throughout the QGP. The "input" charm-quark spectra 
are obtained from either perturbative QCD (pQCD) cal- 
culations or thermal distributions and the inelastic re- 
action employed is the traditional gluo-dissociation pro- 



cess [5]. A narrowing of the J/ip rapidity distribution is 
predicted in A- A relative to p-p collisions. In the statisti- 
cal hadronization model (SHM) [10] all charmonia are pro- 
duced from coalescence of c and c quarks at the hadroniza- 
tion transition. The relative abundance of open and hid- 
den charm states is determined by their mass and spin- 
isospin degeneracy, based on the assumption of thermal 
equilibrium at the hadronization temperature, T c . The 
resulting J jip rapidity distributions from the SHM are 
also narrower in A-A than in p-p collisions, but signifi- 
cantly wider than those from the kinetic recombination 
approach, mostly due to broader input distributions of 
charm-quark cross section. In the comovers interaction 
model (CIM) [11] the primordially produced J ftp's are 
subject to a series of cold nuclear matter induced (initial 
state) effects and subsequently destroyed by "comovers" 
in the medium with an effective dissociation cross section; 
&Diss ~0.65-l mb. The CIM has recently been augmented 
by including charm quark coalescence into J/ip [12]. The 
rapidity dependence of experimental data is reproduced 
both for Au-Au and Cu-Cu collisions at RHIC, mostly 
due to initial-state effects in the incoming nuclei. 

In the present work we apply a rate-equation approach 
5,6^ to 200 AGeV Au-Au collisions at RHIC, to calculate 
J ftp production at forward rapidity and compare to pre- 
viously obtained midrapidity results [6J . The approach as- 
sumes a thermalized medium (QGP and subsequent had- 
ronic gas (HG)) in which anomalous suppression and re- 
generation from coalescence of c-c quarks and D-D mesons 
occur. In Section [2] we evaluate cold nuclear matter ef- 
fects which affect the charmonium abundances before the 
medium thermalizes. We then recall basic ingredients of 
our approach to evaluate J/ip suppression and regenera- 
tion in QGP and HG phase in Section [31 We discuss the 
numerical results for the centrality dependence of the in- 
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elusive J/ip yield and (p 2 ) as well as \c to J/ip and ip' to 
J ftp ratios in Section[4]and the J/ip transverse momentum 
spectra in Section [5] Conclusions are given in Section [6] 



2 Cold nuclear matter effects 

A "pre-charmonium" cc pair produced in an initial hard 
nuclcon-nucleon (N-N) collision first travels through the 
incident (cold) nuclei before becoming a fully formed char- 
monium (W = J/ip, Xci i 3 ')- The modification of initial par- 
ton distribution functions, affecting the hard production of 
pre-charmonium, and the interactions between pre-!? - and 
cold nuclei are referred to as cold nuclear matter (CNM) 
effects. The latter provide a baseline for identifying the 
anomalous suppression and/or enhancement of charmonia 
which is hoped to give insights about the subsequent hot 
medium. An accurate estimate of CNM effects is therefore 
mandatory. In the present work we assume CNM effects to 
be the only relevant ones to charmonium production be- 
fore the medium thermalizes and investigate them within 
the following two schematic baseline scenarios: (1) Nuclear 
absorption + Cronin effect. (2) Shadowing + Nuclear ab- 
sorption + Cronin effect. 

Let us start with scenario 1: The initial dissociation by 
primordial nucleons passing by is described by the stan- 
dard Glauber model resulting in a spatial charmonium 
distribution at the thermalization time To, 

fo(x t ,T ) = (Tp p J d 2 s dz dz'p A {s,z) pEs(x t - s,z') 

x exp j- J dz A p A {s, z A )a n -, 

Xexpj-y dz B pB(x t - S,Z B )(Tnuc\ , (1) 

where pa,s are Woods-Saxon profiles [13] of nuclei A and 
B. The nuclear absorption cross section, cr nuc =1.5 mb 
(oVmc=2.7 mb for tp'), serves as a parameter to regulate 
the strength of dissociation due to all CNM effects com- 
bined, estimated from midrapidity d-Au collisions [14] . 
Shadowing effects are effectively parameterized into a nuc , 
which we furthermore assume to be rapidity-independent. 
Such a scenario may be justified if shadowing in one Au 
nucleus (forward rapidity in d-Au) is roughly compen- 
sated by anti-shadowing in the other Au nucleus (back- 
ward rapidity in d-Au). In addition to absorption, trans- 
verse momentum (p t ) spectra of (prey's are broadened 
compared to p t spectra in p-p collisions due to the Cronin 
effect. We assume the physical mechanism of the Cronin 
effect to be the rescattering of gluons in the CNM be- 
fore they fuse into charmonium. Therefore, the increase 
of (p 2 ) is proportional to the mean length, (I), the gluon 
travels in CNM assuming a random walk treatment of 
the gluon rescattering, resulting in A(pf) = {p 2 ) AA — 
(Pt)pp = a gN(l)- The proportionality coefficient a g N can 
in principle be estimated from d-Au data: here we adopt 
a gN =0A GeV 2 /fm which is compatible with SPS NA50 
data , 15] and current RHIC data [14] . At forward rapidity 



the d-Au data show a pt broadening of A(p 2 )=(p 2 



't/dAu — 

(p 2 ) pp =0.8±0.4 GeV 2 [14] (we combined the uncertainty 
for \p 2 ) of d-Au and p-p in quadrature). With a "path 
length" (Z)=4 fm for d-Au collisions [15] we estimate an 
uncertainty of 0,^=0.1^0.3 GeV 2 /fm in our numerical 
calculations and then we perform a Gaussian smearing 
over the input power-law p t spectra from p-p data with 
A(jft)=a gN (l). 

Concerning scenario 2 we follow the treatment of Ref. 
[12j and separately treat nuclear absorption, shadowing 
and Cronin effect. The total suppression by CNM effects 
(nuclear absorption + shadowing) in this scenario is com- 
parable to the scenario 1 at midrapidity but is stronger at 
forward rapidity due to stronger shadowing at forward ra- 
pidity. Concerning p t distributions and Cronin effect our 
treatment in this scenario is identical to scenario 1. 



3 Kinetic charmonium production approach 

We assume that the medium formed at RHIC reaches ther- 
mal equilibrium at about to=1/3 fm/c and use a thermal 
rate equation thereafter to describe the subsequent evolu- 
tion of 'F's according to, 



~dT 



(2) 



(as usual, we account the feeddown from x c 's and ip"s to 
J ftp's', formation time effects, as considered in Ref. |16| . 
are not included in the present analysis). The initial con- 
dition follows from CNM effects discussed in the previous 
section. The loss term, —r$,N#, accounts for the dissoci- 
ation of the primordially produced charmonium and the 
gain term, r^N^ 1 , for the regeneration of charmonia via 
coalescence of c and c quarks. /V is the in- medium char- 
monia dissociation rate and iV^ q is the equilibrium limit 
of the charmonium abundances. The rate equation treats 
the coalescence of c and c quarks as a continuous process 
over the course of medium evolution in line with lattice 
QCD [T7] and potential-model results [T^ which suggest 
the J/ip to survive at temperatures well above the critical 
temperature T c . 

Let us first address the dissociation of primordially 
produced charmonia (direct component), which is identi- 
fied as the solution of the homogeneous rate equation with 
only the loss term included. In our previously constructed 
kinetic approach a Boltzmann transport equation is em- 
ployed to describe the evolution of the charmonium phase 
space distribution functions f&(x,p,r) in a thermalizcd 
medium with initial distribution from Eq. (fT|). 

P t "d ll fq,(x,p,T) = -E w r#(x,p,T) fy(x,p,r) . (3) 

We neglect elastic charmonium rescattering and only ac- 
count for parton(hadron)-induced inelastic scattering. The 
momentum-dependent dissociation rate is calculated 
via a convolution of the charmonium dissociation cross 
section, a dlss , with a thermal distribution f(u>;T) of par- 
ticles from the heat bath (QGP or HG) and their relative 
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velocity v re i, with charmonia [6J, 

Mx, p, t) = J2J (2^)3 f(^T(r)) 4t Vrel ■ (4) 

In QGP wc include the in-medium reduction of charmo- 
nium binding energy due to color Debye screening. Under 
these circumstances the gluo-dissociation process \P + g — > 
c + c [9] is ineffective to destroy a charmonium [7] and we 
instead employ a quasifree dissociation process: i + W — ► 
i + c + c (i—g,q,q) [7j. The strong coupling constant a s in 
the quasifree cross section is one of the two main param- 
eters in our approach and is adjusted to reproduce the 
J/ip yield in central Pb-Pb collisions at SPS (resulting 
in a s =0.24, cf. the upper panel of Fig. [2]). The compar- 
ison between the charmonium 3-momentum dependence 
of gluo-dissociation (with vacuum binding energy) and 
quasifree rates (with in-medium binding energy) can be 
found in Ref. [BJ: these two dissociation rates are com- 
parable for low-momentum J/tjj but at high momentum 
the gluo-dissociation rate is significantly smaller than the 
quasifree one; for Xc the gluo-dissociation rate is much 
larger than quasifree rate. It is noteworthy that the so- 
called "leakage effect" [TW2"0"Il2l] is accounted for by the 
Boltzmann transport equation approach: high p t charmo- 
nia which travel outside the fireball are not subject to 
suppression anymore rendering less suppression. 

To solve the transport equation (|3|) we need to know 
the medium temperature evolution, T(r), which is deter- 
mined in the following way: We assume the total entropy 
to be produced at the thermalization time (to=1/3 fm/c 
at RHIC). The subsequent expansion is approximated by 
a cylindrically and isentropically expanding fireball, 

V fb {t) = (z + v z t) tt (r + \ a ^ 2 f ■ (5) 

A freeze-out temperature of Tf Q ~ 120 MeV terminates 
the evolution and results in r/ G = 10— 15 fm/c. The fireball 
covers 1.8 rapidity units and the initial transverse radius 
ro represents the initial transverse overlap of the two col- 
liding Au nuclei at a given impact parameter b. The expan- 
sion parameters v z , a z , a ± capture the main aspects of hy- 
drodynamic calculations to reproduce observed flow veloc- 
ities. For simplicity we assume that the expansion param- 
eters are the same for midrapidity (\y\ <0.35) and forward 
rapidity (|y| €[1.2,2.2]). The total entropy is, however, dif- 
ferent, being inferred from number of the charged-particle 
rapidity density, for which we use BRAHMS data [55] 
dN c f l /dy=660 at midrapidity and dN c f l /dy=615 at for- 
ward rapidity (for 0-5% centrality); PHOBOS [23 and 
PHENIX [51] give consistent results with slightly larger 
charged particle multiplicities at both midrapidity and for- 
ward rapidity. For the equation of state (EoS) we use an 
ideal gas of massive quarks and gluons for T > T c and a 
resonance gas equation of state for T <T C including the 37 
lowest lying mesons and 37 lowest lying baryons. The crit- 
ical temperature, T c =180 MeV, at RHIC 200 AGeV is in 
line with thermal-model fits of particle ratios and the cur- 
rent predictions of lattice QCD [55]. Combining entropy 
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Fig. 1. (Color online) Temperature profiles at RHIC 200 
AGeV (solid line: midrapidity; dashed line: forward rapidity) 
and at SPS 17.3 AGeV (dot-dashed line) for central collisions 
(Ap ar t =360) of heavy nuclei. 

conservation and EoS the temperature evolution profile 
follows as a function of time r as displayed in Fig. [T] for 
central Au-Au collisions at both midrapidity and forward 
rapidity, and compared to SPS Pb-Pb ^/s=17.3 AGeV col- 
lisions. We can see that the lifetime of QGP at forward 
rapidity is a slightly shorter than at midrapidity due to 
less produced particles and smaller energy densities. 

Now we turn to the regeneration component which is 
identified as the difference between the solution of the full 
rate equation ([5J and the equation without the gain term, 
r&NfP on r.h.s. The latter is dictated by detailed bal- 
ance with the same (dissociation) rate for the direct 
component. For the charmonium abundances in equilib- 
rium limit, A^ q , we employ the statistical model [TO. with 
N^ q =N^ at TZ, where TZ is a correction coefficient which 
will be detailed later, and A|* at is the abundance from 
the statistical model. This assumes that all the c and c 
quarks are exclusively produced primordially, and popu- 
late open and hidden charm states according to relative 
chemical equilibrium. The equilibrium abundance of re- 
generated charmonia is given by 

N$ tat = -fiVn* (6) 

with V: the volume of the fireball, : chemical equilib- 
rium density of charmonium and 7 C : charm quark fu- 
gacity. The latter is determined by the canonical charm 
conservation equation, 

at 1 t/ h(lcVn p) 2 

Ace = olcVn op ——— + "f c Vn hld 7 

with N c5 : total number charm quark pairs produced, n op : 
density of all open charm states, rihid- density of all hid- 
den charm states. In the present work we use the charm 
quark cross section da ec /<iy(y — 0)=95 /ib in line with 
recent PHENIX measurements [UJ (with the assumption 
that 30% of the total produced charm are uniformly dis- 
tributed over 1.8 rapidity units which is the rapidity cov- 
erage of our fireball) . For the rapidity dependence we take 
guidance from perturbative QCD calculation [57], leading 
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to d<r gc /dy(y=1.7)=60 ph. Currently the uncertainties of 
charm production in both theoretical predictions and ex- 
perimental measurements are rather large. The modified 



Bessel function factor 



hhcVn op ) 



on the r.h.s accounts for 



IoilcVriap) 

exact cc conservation in the canonical ensemble relative 
to grand-canonical ensemble [2"5lf2"9"] . 

The abundance of charmonia, N^ tat , from the statisti- 
cal model is subject to two additional corrections to obtain 
in the rate equation Eq. © . The first correction sim- 
ulates incomplete charm quark thermalization in medium: 
It is natural to expect that the coalescence rate from 
non-thermalized c and c quarks is smaller than for fully 
thermalized charm quarks [30il3l] . We implement this cor- 
rection by multiplying the charmonium abundances from 
the statistical model in a schematic way by a factor 1Z = 
1 — cxp(— t/t^j), where r^ q is the thermal relaxation time 
of charm quarks which is the second main parameter in our 
approach. The second correction is a "correlation volume" 
V corr for the coalescing charm quarks [5]. The correlation 
volume increases W production because the cc pairs do not 
have time to populate the entire fireball volume. We im- 
plement this correction by replacing the fireball volume 
V in the argument of the Bessel function by a correla- 
tion volume V corr in Eq.([7|). The "correlation volume" is 
estimated as the volume explored by a receding cc pair: 
Vo(T)=47r(r + (w c )r) 3 /3, where ro— 1.2 fm represents a 
minimal radius characterizing the range of strong interac- 
tions, and (v c ) denotes the average relative speed of the 
produced c and c quarks which we vary between 0.5c and 
0.8c. Within this setup we adjust our second parameter 

to the inclusive J ftp yield in central Au-Au at RHIC 
at midrapidity, resulting in Tg (J =4-7 fm/c, corresponding 
to (w c )=0.8c and 0.5c, respectively. 

The transverse momentum distributions of regener- 
ated charmonia are approximated by local thermal equi- 
librium and boosted by the transverse flow of the medium, 
corresponding to a blastwave expression 32J, 



Ptdpt 



ex rat 



coal 



rdrKi 



( m t cosh y t 



T 



p t sinh y t 



T 



with mt=y/rn# + pf, transverse rapidity y t =tanli 

r 

' }s R 



(8) 
_1 v t (r) 

with linear flow profile v t (r)=v s j^ and surface velocity v s 
given by the fireball evolution formula, Eq.©. We evalu- 
ate the blastwave expression at the hadronization transi- 
tion (T c =180 MeV) and neglect rescattering of tf^'s in the 
hadronic phase. 



4 Centrality and Rapidity Dependence 

The number of J/ip produced in initial hard collisions at 
given impact parameter (6) is estimated from the pro- 
duction cross section in p-p collisions and then scaled by 
number of binary collisions in N-N collisions, N co u(b). 
We use the J/ijj production cross section from PHENIX 

data 33 with d<Tp/ ,/ '/d?/=750(500) nb (with about 30% 
uncertainty) at mid (forward) rapidity. The centrality de- 
pendence of the inclusive J/tj) yield in terms of nuclear 
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approach for Rj/a vs 



(Color online) Results of the thermal rate-equation 
centrality at SPS (upper panel) and 
RHIC midrapidity (lower panel) are compared to NA50 [2] and 
PHENIX data ;3 . Solid line: total J/tp yield; dashed line: sup- 
pressed primordial production; dot-dashed line: regeneration 
component; dotted line: primordial production with nuclear 
absorption only. 



modification factor, 



Raa = 



N AA 



N P j%N coll 



(9) 



is displayed in Fig. [2j For central collisions the direct and 
regeneration component are about equal (quite similar to 
Ref. [31]), while for peripheral collisions the direct com- 
ponent dominates because the lifetime of the medium is 
shorter so that there is less time for regeneration pro- 
cesses and for charm quarks to thermally equilibrate (re- 
ducing coalescence into charmonia). Fig. [5] also includes 
the strength of suppression caused by CNM effects (dot- 
ted line) and by the hot medium indicating the latter to 
be substantially larger than the former. 

Next we proceed to the inclusive J / ip yield at forward 
rapidity as shown in Fig. [3] with scenarios 1 and 2 for CNM 
effects. In scenario 1, where the CNM induced suppression 
at forward rapidity is assumed to be the same as at mid- 
rapidity, the hot medium causes slightly less suppression 
than at mid-rapidity due to a shorter QGP lifetime (recall 
Fig. [1]). Concerning the regeneration component, we note 
that the number of charm quark pairs at RHIC energy is 
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Fig. 3. (Color online) Results of the thermal rate-equation ap- 
proach for Rjjf vs. centrality at forward rapidity compared to 
PHENIX data 3! with CNM effects in scenario 1 (upper panel) 
and scenario 2 (lower panel). Solid line: total J/ip yield; dashed 
line: suppressed primordial production; dot-dashed line: regen- 
eration component; dotted line: primordial production with 
CNM effects only. 



< i 
< 



| 0.6 
< 
< 

OS 0.4 



0.2 



Au-Au (200 A GeV) □ phenk 

— w/ shadowing 

— w/o shadowing 



50 



100 



150 200 

N , 
part 



250 300 350 400 



Fig. 4. (Color online) The thermal rate-equation approach 
results of ratio of forward and midrapidity J/lp Raa versus 
N par t with cold nuclear matter effects in scenario 1 (solid line) 
and scenario 2 (dashed line) compared to PHENIX data [3]- 



between the canonical and grandcanonical limit, so that 
]\[stat ^ jya_ w j^] 1 a between 1 (canonical limit) and 2 
(grandcanonical limit) as following from the charm con- 
servation equation |(7J). If charm and J/ip production in p-p 
collisions (the denominator in Raa) are less at forward ra- 
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Fig. 5. (Color online) Results of the thermal rate-equation 
approach for the (xd + Xc2)/{J/ip) ratio (upper panel) and 
ip'/(J/ip) ratio (lower panel) vs. centrality at midrapidity at 
RHIC. Solid line: the results of the thermal rate-equation ap- 
proach; dashed line: the nuclear absorption only. The J/ip's in 
the denominator include feeddown from \ c and Va- 



pidity than at midrapidity by about the same fraction, the 
reduction of regeneration component at forward rapidity 
is about Raa ~ A^ -1 . Adding up the two components the 
total J/ip yield at forward rapidity is almost equal to that 
at midrapidity since the slight decrease in the regenera- 
tion component is compensated by the slight increase in 
the direct component, see upper panel of Fig. [3] Scenario 
2 differs from scenario 1 in that the direct component at 
forward rapidity is subject to stronger CNM suppression. 
Consequently, the inclusive J / ip yield at forward rapidity 
is more suppressed which leads to better agreement with 
experimental data as seen in the lower panel of Fig. [3] 
To better illustrate the comparison between forward and 
midrapidity we plot the ratio of forward and midrapidity 
J/ip Raa and compare to experimental data in Fig. [4] 
Even with a stronger cold matter induced suppression at 
forward rapidity due to shadowing (scenario 2) our ap- 
proach still appears not to fully account for the experi- 
mental findings. However one needs to keep in mind that 
the current uncertainties in several key inputs, in partic- 
ular charm production cross section at forward rapidity 
and shadowing, are still appreciable. More accurate d-Au 
data will enable more definite conclusions. 

Finally we calculate the Xc/iJ/ip) and ijj' / (J/ip) ratios 
which may provide additional discrimination of charmo- 
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nium production mechanisms: e.g., the gluo-dissociation 
process and quasifree dissociation process give compara- 
ble suppression for J/tp but for \c the gluo-dissociation 
process gives much larger dissociation rate [6] leading to a 
much smaller Xc/(J/ip) ratio. As another example, if for- 
mation time effects [20ll34ll35] are important one may ob- 
serve less suppression of Xc than J jip> due to a longer 35J 
formation time of Xc compared to J /if), together with a 
smaller dissociation cross section for a "pre-hadronic" cc 
pair than for a fully formed charmonium. This is in con- 
trast to most standard dissociation mechanisms (since Xc 
is a higher excited cc state than J /if) with smaller binding 
energy and therefore is more easily destroyed) and regen- 
eration mechanism (x c IS heavier so that its equilibrium 
abundance is suppressed compared to J /if) by the Boltz- 
mann thermal factor). 

For Xc states, we constrain ourselves to Xci and Xc2 
with a combined average branching ratio of 27% into J ftp's 
(36] . The resulting Xc/ {J/i>) an d fp' / (J/ip) ratio at midra- 
pidity are shown in Fig. [5] Both drop with centrality be- 
low the ratios obtained from CNM-induced suppression 
(we assume that J ftp and Xc undergo the same strength 
of CNM-induced suppression). 



5 Transverse Momentum Spectra 

In order to gain further insights into charmonium produc- 
tion mechanisms, and to discriminate among them, we 
investigate the transverse momentum spectra of J/ip. Let 
us again start at midrapidity: The centrality dependence 
of (pf ) is shown in the upper panel of Fig. [6j The (pf ) for 
the direct component increases with centrality due to a 
stronger Cronin effect generated by a longer path length 
(I) travelled by the gluon pair (recall Section [2|) in more 
central collisions. The (pf ) of the regeneration component 
also increases with centrality due to an increase of the ra- 
dial flow velocity at the end of mixed phase. However, the 
magnitude of (pf ) of the regeneration component is much 
smaller than primordial production at all centralities. The 
main point is that for more central collisions the regener- 
ation component makes up an increasing fraction of the 
total yield leading to an almost constant (pf ) with central- 
ity, consistent with current PHENIX data. Similar results 
are also obtained by the kinetic recombination model of 
Ref. [8j. On the contrary, at SPS energy, the direct compo- 
nent dominates at all centrality so that the Cronin effect 
dominates leading to a monotonically increasing {pi), see 
lower panel of Fig. [6j A decreasing (p|) with centrality 
is thus a supporting signature of the presence of regen- 
eration. More definite conclusions will be possible with 
improved data accuracy. 

We now turn to the (pf) vs. centrality at forward ra- 
pidity as shown in Fig. [7] (with CNM effects in scenario 
1). Within the current uncertainties in the Cronin effect 
our results agree with experimental data. The increase of 
(pf ) suggests that the direct component at forward rapid- 
ity is substantial; otherwise the (pf ) would be significantly 
smaller. 
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Fig. 6. (Color online) Results of the thermal rate-equation 
approach for (pf) vs. centrality at RHIC midrapidity (upper 
panel) and SPS (lower panel) compared to PHENIX [5] and 
NA50 data |37lll5] . Solid line: (pf) for the total J/ip yield (di- 
rect component + regeneration component); dashed line: (pf) 
for the direct component; dot-dashed line: (pf) for the regen- 
eration component. 
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Fig. 7. (Color online) Results of the thermal rate-equation 
approach for (pf) vs. centrality at forward rapidity compared 
to PHENIX data [3]. The band represents the (pf) for the to- 
tal J ftp yield (direct + regeneration). The width of the band 
represents the uncertainty from the strength of the Cronin ef- 
fect (estimated from d-Au data with a g N= 0.1-0.3 GeV 2 /fm). 
The cold nuclear matter effects are implemented in scenario 1. 
Dot-dashed line: (pf ) for the regeneration component. 
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Fig. 8. (Color online) Results of the thermal rate-equation approach for R'aa vs. transverse momentum for different centrality 
selections at midrapidity are compared to PHENIX data 3 . Solid line: total J/ip yield; dashed line: suppressed primordial 
production; dot-dashed line: regeneration component. The Cronin effects are implemented with a g N= 0.1 GeV 2 /fm. 
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Fig. 9. (Color online) Results of the thermal rate-equation approach for R^a vs - transverse momentum for different centrality 
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production; dot-dashed line: regeneration component. The CNM effects are implemented in scenario 1 with a g N= 0.2 GeV 2 /fm. 
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For a more differential comparison with experimental 
data we plot R-AAipt) for different centrality selections at 
both mid- and forward rapidity (with CNM effects in sce- 
nario 1 and a g jv=0.2 GeV 2 /fm) in Fig. [8] and [9] respec- 
tively. Overall, our model results reproduce B,AA(pt) data 
reasonably well. By comparing forward and midrapidity 
results there is an indication for stronger suppression at 
forward rapidity relative to midrapidity mainly in the low 
Pt regime (for central and scmiccntral collisions). Less re- 
generation and stronger Cronin effect at forward rapidity 
can both generate such an effect. 



6 Conclusion 

In this work we have applied a previously constructed 
thermal rate-equation approach to study charmonium pro- 
duction in 200 AGeV Au-Au collisions at RHIC at forward 
rapidity and compared to midrapidity results. Our calcu- 
lations reasonably reproduce the experimental data of in- 
clusive J ftp Raa, (j>t) an d Pt spectra at forward rapidity. 
We find that the observed stronger suppression at forward 
rapidity can be partially attributed to smaller regenera- 
tion resulting from smaller open charm production but 
additional suppression from cold nuclear matter effects 
(shadowing) is required as well. An accurate knowledge 
of shadowing will be needed to clarify the "puzzle" of the 
stronger suppression at forward rapidity. As a next step 
we will improve our current treatment of the coalescence 
of c and c quarks by performing a microscopic calculation 
with a time-dependent charm quark distribution based on 
Ref. [38] . which will allow for a more accurate evaluation 
of regeneration effects, especially their dependence on in- 
complete charm-quark equilibration. We also look forward 
to new experimental data to better constrain the inputs 
of our approach which will pave the way for more quanti- 
tative conclusions. 
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